ON THE PERIODIC SOLUTIONS 
OF A PERTURBED DOUBLE PENDULUM 

JAUME LLIBREi, DOUGLAS D. NOVAES^ AND MARCO ANTONIO TEIXEIRA^ 
To Waldyr Oliva in his 8(f^ Birthday 



Abstract. We provide sufficient conditions for the existence of periodic solu- 
tions of the planar perturbed double pendulum with small oscillations having 
equations of motion 

e'l = -2a6»i +06*2 +£Fi(t, 91,^1, 6*2, 62), 

6»2 = 2ae-i -206*2 +eF2{t, 6*1, ei, 6*2, 02), 
where a and e are real parameters. The two masses of the unperturbed double 
pendulum are equal, and its two stems have the same length I. In fact a = g/l 
where g is the acceleration of the gravity. Here the parameter e is small and 
the smooth functions Fi and F2 define the perturbation which are periodic 
functions in t and in resonance p:q with some of the periodic solutions of the 
unperturbed double pendulum, being p and q positive integers relatively prime. 



1. Introduction and statement of the main results 

We consider a system of two point masses mi and TO2 moving in a fixed plane, in 
which the distance between a point (cahed pivot) and mi and the distance between 
mi and m2 are fixed, and equal to li and I2 respectively. We assume the masses 
do not interact. We allow gravity to act on the masses mi and m2. This system is 
called the planar double pendulum. 

The position of the double pendulum is determined by the two angles 9i and 62 
shown in Figure 1. We consider only the motion in the vicinity of the equilibrium 
6*1 = ^2 = 0, i.e. we are only interested in small oscillations around this equilibrium. 
Expanding the Lagrangian of this system to second order in 9i and 62 and their 
time derivatives, the corresponding Lagrange equations of motion are 

{mi + m2)li9i + 7712^2^2 + {mi + TO2).9^i = 0, 

TO2^2^2 + m2ll9i + 77125^2 = 0. 

For more details on these equations of motion see [2] . Here the dot denotes deriv- 
ative with respect to the time t. 

If we take mi = TO2 — m > 0, li = I2 = I > Q and divide the above equations by 
ml, system ([T]) becomes 

2 ^'1 + ^2 + 2^01 =0, 

(2) .. .. ; 

01+02 + 7^2 =0. 
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Figure 1. The planar double pendulum. 

Isolating 9i and 62 of equations ([2]) and denoting a = g/l, we obtain that the 
equations of motion of the double pendulum with small oscillations, two equal 
masses m and having the two stems the same length / are 

01 = -2a0i+a92, 

(3) 

02 = 206*1 - 2a02, 
where g is the acceleration of the gravity. 

Taking a new time t given by the rescaling r — ^/a t the equations of motion Q 
become 

0'i^-20i+02, 

(4) 

02 — 201 — 202, 

where now the prime denotes derivative with respect to the new time r. 

The objective of this paper is to provide a system of nonlinear equations whose 
simple zeros provide periodic solutions of the perturbed planar double pendulum 
with equations of motion 

0'; = -20i + 02+eFi{T,0i,0[,02,0'2), 

(5) 

02 = 201 - 202 + eF2{T,9u9[,02,9'2), 

where e is a small parameter. Here the smooth functions Fi and F2 define the 
perturbation. These functions are periodic in r and in resonance p:q with some of 
the periodic solutions of the unperturbed double pendulum, being p and q positive 
integers relatively primes. In order to present our results we need some preliminary 
definitions and notations. 

The imperturbed system ^ has a unique singular point, the origin with eigen- 
values 



±y2~V2i, ±y2 + V2i. 

Consequently this system in the phase space {0i,0[,02,02) has two planes filled of 
periodic solutions with the exception of the origin. These periodic solutions have 
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periods 



27r 27r 
Ti = or T2 = , 

V2-\/2 V2W2 



according they belong to the plane associated to the eigenvectors with eigenvalues 
±\/2 — \/2i or ±^72 + i, respectively. We shall study which of these periodic 
solutions persist for the perturbed system ([5]) when the parameter e is sufficiently 
small and the perturbed functions Fi for i = 1, 2 have period either pTi/ q, or PT2/ q, 
where p and q are positive integers relatively prime. 

We define the functions: 



j^^ ' sin (^V^2- V2t^ (V2F1 + F2) dr, 



ai(Xo,ro) 

(6) 

^2(^0, Fo) - ' COS (^\/2- V2r^ (72 Fi + F2) dr, 



where 

(7) Ffc = Ffc(T,Ai,Bi,Ci,7?i) 

for fc = 1 , 2 with 

Ai = , fXo cos (72^^72 +yosinfv/2-\/2r 

v/4- 2V2 ^ ^ ^ ^ 

Bi = (yoCos(V2-\/2r) - Xo sin ( ^2 - \/2 r) ) , 

Ci = , (Xo cos fv/2-\/2r') + Yo sin (^2-^2^ 

v/2- V2 ^ ^ ^ ^ 



Di = cos ( V2W2 r — Xo sin ( v/2- V2 r 
A zero (Xq , ^o*) '^^ t'^c nonlinear system 

(8) ei(Xo,ro) = o, e2(^o,io) = 

such that 

'9(^1,^2) 



det 



^0, 



(Xo,Yo) = (X-,Yo-), 



^9(Xo,yo) 

is called a simple zero of system ([8]). 

Our main result on the periodic solutions of the perturbed double pendulum ([5]) 
which bifurcate from the periodic solutions of the unperturbed double pendulum 
(HI) with period Ti traveled p times is the following. 

Theorem 1. Assume that the functions Fk of the perturbed double pendulum with 
equations of motion ^ are periodic in r of period pTi/q with p and q positive 
integers relatively prime. Then for s sufficiently small and for every simple 
zero (Xo,yo*) 7^ (Oi 0) of the nonlinear system (|8]), the perturbed double pendulum 
([5|) has a periodic solution (6'i(r, e), 6'2(t, e)) tending when e to the periodic 
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solution {9i{t),62{t)) run p times of the unperturbed double pendulum (j4]) given by 
^ (^X* cos (V2-\/2r) + Yq* sin (^2- \/2i 



(9) 



V4 - 2V2 



:!= (x* COS (^2- V2r) + Fq* sin (^2- V2r))^ 



Theorem [T] is proved in section [2] Its proof is based in the averaging theory for 
computing periodic solutions, see the appendix. 

We provide an application of Theorem [T] in the following corollary, which will be 
proved in section |31 



Corollary 2. //F2(6'i, 6li, 6*2, 6*^) = (l - 6*? ) sin (^^2 - V^tj and i^i(6'i, 6''i, 6*2, 6*^) = 

0, then the differential equation ([5]) for e ^ sufficiently small has two peri- 
odic solutions {9i(T,e),92{T,e)) tending when e — )■ to the two periodic solutions 
{9i(t),92{t)) of the unperturbed double pendulum Q given by ^ with 



(X;,Y;)= 2,/2 ,0 and (X„-,r„*) = 



respectively. 

Now we define the functions 

rpT: 



(10) 



g\Zo,WQ) = £ sin (^^J2 + V2t^ (-\/2 Fi + -^2) dr, 



g^iZo, Wo) = £ ' COS (^V^2 + V2r^ (-V2F1 + F2) dr, 
where Ffe = ^^(t, A2, B2, C2, 1)2) for fc = 1, 2 with 

A2 = , ^ fZoCosfV2 + \/2T) + W^osinfV2 + \/2r)) 

V4 + 2V2^ V y V yy 

= (-M^ocos(V2T72r) + Zo sin ( ^2 + %/2 t) ) , 

C2 = , f Zo COS (V2 + V2t) + Wo sin (^2 + V^r 

V2 + V2 ^ ^ ^ ^ 



D2 = W^o COS (^V2 + V2 tJ - Zo sin (^V2 + V2r 
Consider the nonlinear system 

(11) g\Zo,Wo) = o, g^{Zo,Wo) = o. 

Our main result on the periodic solutions of the perturbed double pendulum ([S]) 
which bifurcate from the periodic solutions of the unperturbed double pendulum 
Q with period T2 traveled p times is the following. 

Theorem 3. Assume that the functions Fk of the perturbed double pendulum with 
equations of motion ^ are periodic in r of period pT2/q with p and q positive 
integers relatively prime. Then for e ^ sufficiently small and for every simple 
zero {Zq, Wq) ^ (0,0) of the nonlinear system ([TT|) . the perturbed double pendulum 
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([5|) has a periodic solution (9i{T,e),02{T,e)) tending when e — > to the periodic 
solution {0i{t),62{t)) run p times of the unperturbed double pendulum (j4]) given by 




Theorem |3] is also proved in section [51 

Now we provide an application of Theorem |3] in the next corollary, which will be 
proved in section |31 

Corollary 4. //i^i(6'i, 6li, 6*2, 6*^) = 0'2 + ef cos (^\/ 2 + V2 1^ and F2iei,9[,e2,9'2) ^ 
0, then the differential equation ([5]) for e 7^ sufficiently small has one periodic so- 
lution (6'i(t, e),6'2(T, e)) tending whene to the periodic solution (6'i(t), 6'2(t)) of 
the unperturbed double pendulum ^ given by (da) with {Z^, W^) = (0, -8 (2 + V2)) . 



2. Proof of Theorems [T] and [3] 

Introducing the variables {x,y,z,w) = {61, 9'^, 92 ,62) we write the differential 
system of the perturbed double pendulum ([S]) as a first-order differential system 
defined in R^. Thus we have the differential system 

—2x + z + eFi(r, a:, y, w), 
2x — 2z + eF2{T, X, y, z, w). 



(13) 



y 

z' 
w' 



System (I13|) with e = is equivalent to the unperturbed double pendulum system 
^ , called in what follows simply by the unperturbed system, otherwise we have the 
perturbed system. 

We shall write system (|13p in such a way that the linear part at the origin will be 
in its real normal Jordan form. Then, doing the change of variables (r, x, y, z, w) — >■ 
(T,X,y,Z, W) given by 
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the differential system (|T3|) becomes 

X' = V2-\/2y, 

Y' = -v/2- + e]- (^/2 Fi + F2 
(15) ^ 

where f:,(T, X, Y, Z, W) = F,{t, A, B, C, 2?) for i = 1,2 with 



C = 



V/4-2V2 ^2 (2 + V2) ' 



V2-\/2 V2 + \/2' 
2? = r + VF. 

Note that the hnear part of the differential system ()15p at the origin is in its real 
normal Jordan form. 



Lemma 5. The periodic solutions of the differential system (jlSp with e = are 



X{t) ==XoCos(V2-\/2t) +rosin(V2-\/2T 



(16) F(r) =roCos(v/2- V2rj -Xosin(v/2- V2Tj , 

Z(t) = 0, 
W{t) = 0, 

0/ period Ti , and 

^(r) -0, 

y{r) = 0, 

(!'') Z{t) = Zocos (V2 + \/2r) + Wo sin {\f2^W2T 



W{t) = Wq cos V2+72 rj - sin [^2 + ^/2 rj , 

0/ period T2 . 

Proof. Since system p5|) with e = is a linear differential system, the proof follows 
easily. □ 



Proof of Theorem QJ Assume that the functions Fk of the perturbed double pen- 
dulum with equations of motion ^ are periodic in r of period pTi/q with p and 
q positive integers relatively prime. Then we can think that system ([5]) is periodic 
in T of period pTi . Thinking in this way the differential system and the periodic 
solutions (fT6|) have the same period pTi . 
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We shall apply Theorem |6] of the appendix to the differential system ([T5|) . We 
note that system (|15|) can be written as system ([T9l) taking 

V2T72W, 



/ ^ \ 

Y 
Z 

T/f/ / 



t = T, Go(t,x) = 



Gi(t,x) = 



1 



and 6*2 (i, x, e) = 



( \ 




V ) 



We shall study which periodic solutions p6)) of the unperturbed system p5)) 
with e = can be continued to periodic solutions of the unperturbed system (jlSp 
for e 7^ sufBciently small. 

We shall describe the different elements which appear in the statement of The- 
orem |6] in the particular case of the differential system (fT5|) . Thus we have that 
17 — R**, k — 2 and n = A. Let ri > be arbitrarily small and let r2 > be 
arbitrarily large. We take the open and bounded subset V of the plane Z — W — 



V = {(Xo, Yo, 0, 0) e R4 : ri < ^X^ + Y^^ < r^}. 

As usual C\{V) denotes the closure of V. If a = {Xo,Yo), then we can identify V 
with the set 



here 



{a e : ri < ||a|| < r2}, 
denotes the Euclidean norm of R^. The function /3 : Cl{V) R^ 



IS 



/3(a) — (0, 0). Therefore, in our case the set 

Z = {z, - {a, /3(a)) , a € C\{V)} = {(^0,^0, 0, 0) G R^ : ri < ^Jx^ + Y^ < r^}. 

Clearly for each Zq S Z we can consider the periodic solution x(r, Zq.) = {X{t),Y{t), 
0, 0) given by (fTB)) of period pTi. 



Computing the fundamental matrix M^,^ (r) of the linear differential system (jlSp 
with £ = associated to the T-periodic solution Zq = {Xo,Yo,Q,0) such that 
A^zc (0) be the identity of R", we get that M (r) = (t) is equal to 



f cos (^^2 - V2t^ sin(^y2~~72r 
-sin(^V2~~72 



2- V2rj 

cos(^V2 + \/2rj sin (^^2 + \/2r^ 

-sin fV2 + \/2t') cosfV2 + \/2r 
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Note that the matrix M^^ (r) does not depend of the particular periodic solution 
x(r, Zq). Since the matrix 



/ Q \ 



2sin^(y2 7r) sin (2^2 tt) 

VO -sin(2V27r) 2sin^ (\/2 7r) / 



satisfies the assumptions of statement (ii) of Theorem [5] because the determinant 

: 4sin2 (\/27r) ^ 0, 



2sin^(V2 7r) sin(2V2 7r) 
-sin(2\/2 7r) 2sin^(\/2 7r) 



we can apply this theorem to system (|15p . 

Now ^ : R4 R2 is ^(X, y, Z, W) = {X, Y). We calculate the function 

/ 1 rpTi ^ 
pTi Jo 



and we obtain 



^C„'MGi(T,x(T,z„))dr , 



G2{Xo,Yo) 



47r 



V2t\ (^2 



F1+F2] dr 



V 



47r 



2-V2t] (V2Fi+F2]dT 



J 



where the functions of Fk for fc = 1, 2 are the ones given in ([Tj- Then, by Theorem 
ini we have that for every simple zero [X^^Y^] S V of the system of nonlinear 
functions 

(18) gi(Xo,ro) = 0, g2[Xo,Ya) = Q, 

we have a periodic solution {X, Y, Z, W){t, e) of system such that 

{X, Y, Z, W){t, e) {X*,Y*,0, 0) as e ^ 0. 

Note that system (ITS)) is equivalent to system because both equations only 
differs in a non-zero multiplicative constant. 

Going back through the change of coordinates (fT4l) we get a periodic solution 
(x, y, z, 'w){t, e) of system (jl6p such that 



/ a;(T,e) \ 
y{T,e) 
z{T,e) 



as e — 0. 



^ (^0 cos V2r) + Y^ sin (^2 - %/2t) ) ^ 

(Yq* cos (V2-\/2r) - X^ sin (72^72 r)) 
cos (V2- %/2t) + Fq* sin (^2- %/2t 



V2-\/2 

Fo* cos (^2- V2r) - X^ sin (^2- V2r) 
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Consequently we obtain a periodic solution 6'2)(t, e) of system ^ such that 



I ^ [X*^ cos (V2-\/2r) + sin (^2 - \/2t) ) \ 



V 



[x^ cos (V2^V2 + Fq* sin (\/2-\/2: 



V2-\/2 

as e — > 0. Hence Theorem [1] is proved. □ 



Proof of Theorem O This proof is completely analogous to the proof of Theorem 
[H □ 

3. Proof of the two corollaries 

Proof of Corollary [H Under the assumptions of Corollary [2] the nonlinear system 
© becomes 

7r(x2 + 3if + 8(-2 + V2)) 



gi(Xo,yo) = — 



^2(^0,^0) 



8 (2 - %/2) 



3/2 



4(2-V2)'/'' 



This system has the following four solutions 

{X*.Y*) = (^2^2 (2-^/2),0^ and (Xo*,^;) = |^0, ±2^1(2^7^) j . 

But the solutions which differs in a sign are different initial conditions of the same 
periodic solution of the unperturbed double pendulum. Moreover, it easy to check 
that these solutions are simple. So, by Theorem [1] we only have two periodic 
solutions of the perturbed double pendulum. This completes the proof of the corol- 
lary. □ 

Proof of Corollary Under the assumptions of Corollary |4] the nonlinear system 
([Tot becomes 



^ (^^^(10-7^2) (2 + V2)Wo - 8^ Zo 



4J2 (2 + V2) 



r2iy w^ ^ {^Wj - 2Wi - 161^0 + iV2Zi - QZj) 
y [Zq, Wq) = 



8y'2 (2 + ^/2) 

This system only has the following two real solutions 

(iZ*,W^-) - (0, -8 (2 + ^/2)) and {Z*,W^) = (0, 0). 

the other two solutions are not real. But the solution (0,0) is not valid. Therefore, 
by Theorem|3]we only have one periodic solution of the perturbed double pendulum. 
This completes the proof of the corollary. □ 
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Appendix: Basic results on averaging theory 

In this appendix we present the basic resuh from the averaging theory that we 
shaU need for proving the main results of this paper. 

We consider the problem of the bifurcation of T-periodic solutions from differ- 
ential systems of the form 

(19) k{t) = Go(t,x) +eGi(t,x) +e2G2(t,x,e), 

with e = to e 7^ sufficiently small. Here the functions Gq , Gi : M x — > M" and 
G2 : M X X (— £o,£o) I^" are functions, T-periodic in the first variable, and 
is an open subset of M". The main assumption is that the unperturbed system 

(20) x(i) = Go(i,x), 

has a submanifold of periodic solutions. A solution of this problem is given using 
the averaging theory. 

Let x(i, z, e) be the solution of the system ([20]) such that x(0, z, e) ~ z. We write 
the linearization of the unperturbed system along a periodic solution x(t, z, 0) as 

(21) y = i?xGo(t,x(i,z,0))y. 

In what follows we denote by Mz (t) some fundamental matrix of the linear differ- 
ential system and by ^ : R*^ x ]R"-'= R'' the projection of M" onto its first 
k coordinates; i.e. ^(xi, . . . , x„) = (xi, . . . ,Xk). 

We assume that there exists a fc-dimensional submanifold Z oi fl filled with 
T-periodic solutions of (|20|) . Then an answer to the problem of bifurcation of 
T-periodic solutions from the periodic solutions contained in Z for system (|19p is 
given in the following result. 

Theorem 6. Let V be an open and bounded subset o/R*^, and let (3 : G\{V) — !• R"~'^ 
be a function. We assume that 

(i) Z ~ {za = (a, P{ol)) , a G Cl(l^)} C and that for each z^ G Z the solu- 
tion x.{t,Za) of ()20|) is T-periodic; 

(ii) for each Za € Z there is a fundamental matrix M^^ (t) of (|2ip such that 
the matrix (0) — (T) has in the upper right corner the k x (n — k) 
zero matrix, and in the lower right corner a {n — k) x {n ^ k) matrix Aq, 
with det(Aa) ^ 0. 

We consider the function Q : Cl(V^) -> R*^ 

(22) g(a)=S,(^j^ M~^Ht)Gi{t,Mt,2a))dt 

If there exists a £ V with g{a) = and dct {(dQ / da) (a)) j^O, then there is a 
T-periodic solution ip(t,e) of system (|19p such that ip{0,e) — > Zq as e — )• 0. 

Theorem IH] goes back to Malkin [3] and Roseau [3] , for a shorter proof see [T] . 
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